arXiv:chao-dyn/9808011v3 14 Sep 1999 


Anomalous scaling regimes of a passive scalar advected by the 

synthetic velocity field 


N. V. Antonov 

Department of Theoretical Physics, St Petersburg University, Uljanovskaja 1, 
St Petersburg, Petrodvorez, 198904 Russia 


Abstract 


The field theoretic renormalization group (RG) is applied to the problem of 
a passive scalar advected by the Gaussian self-similar velocity field with finite 
correlation time and in the presence of an imposed linear mean gradient. 
The energy spectrum in the inertial range has the form E(k ) oc k 1 ~ £ , and 
the correlation time at the wavenumber k scales as k~ 2+v . It is shown that, 
depending on the values of the exponents e and rj, the model in the inertial- 
convective range exhibits various types of scaling regimes associated with the 
infrared stable fixed points of the RG equations: diffusive-type regimes for 
which the advection can be treated within ordinary perturbation theory, and 
three nontrivial convection-type regimes for which the correlation functions 
exhibit anomalous scaling behavior. Explicit asymptotic expressions for the 
structure functions and other correlation functions are obtained; they are 
represented by superpositions of power laws with nonuniversal amplitudes and 
universal (independent of the anisotropy) anomalous exponents, calculated to 
the first order in e and g in any space dimension. These anomalous exponents 
are determined by the critical dimensions of tensor composite operators built 
of the scalar gradients, and exhibit a kind of hierarchy related to the degree of 
anisotropy: the less is the rank, the less is the dimension and, consequently, 
the more important is the contribution to the inertial-range behaviour. The 
leading terms of the even (odd) structure functions are given by the scalar 
(vector) operators. For the first nontrivial regime the anomalous exponents 
are the same as in the rapid-change version of the model; for the second they 
are the same as in the model with time-independent (frozen) velocity field. 
In these regimes, the anomalous exponents are universal in the sense that 
they depend only on the exponents entering into the velocity correlator. For 
the last regime the exponents are nonuniversal (they can depend also on the 
amplitudes); however, the nonuniversality can reveal itself only in the second 
order of the RG expansion. A brief discussion of the passive advection in 
the non-Gaussian velocity field governed by the nonlinear stochastic Navier- 
Stokes equation is also given. St Petersburg University Preprint SPbU IP-98- 
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I. INTRODUCTION 


The investigation of intermittency and anomalous scaling in fully developed turbulence 
remains one of the major theoretical problems. Both the natural and numerical experiments 
suggest that the deviation from the predictions of the classical Kolmogorov-Obukhov theory 
is even more strongly pronounced for a passively advected scalar held than for the velocity 
held itself; see, e.g., JT|-^| and literature cited therein. At the same time, the problem 
of passive advection appears to be easier tractable theoretically: even simplified models 
describing the advection by a “synthetic” velocity held with prescribed Gaussian statistics 
reproduce many of the anomalous features of genuine turbulent heat or mass transport 
observed in experiments, see p|- PH- Therefore, the problem of a passive scalar advection, 


being of practical importance in itself, may also be viewed as a starting point in studying 
anomalous scaling in the turbulence on the whole. 

Recently, a great deal of attention has been drawn by a simple model of the passive scalar 
advection by a self-similar Gaussian white-in-time velocity held, the so-called “rapid-change 
model,” introduced by Kraichnan |T(|; see [8-30] and references therein. For the hrst time, 
the anomalous exponents have been calculated on the basis of a microscopic model and 
within regular expansions in formal small parameters. Within the “zero-mode approach” to 
the rapid-change model, developed in 00 , nontrivial anomalous exponents are related to 
the zero modes (homogeneous solutions) of the closed exact equations satished by the equal¬ 
time correlations. In this sense, the model is “exactly solvable.” The anomalous exponents 
are universal, i.e., they depend only on the space dimension and the exponent entering into 
the velocity correlator. 

Of course, the Gaussian character, isotropy, and time decorrelation are strong departures 
from the statistical properties of genuine turbulence. One step toward the construction of a 
more realistic model of passive advection is the account of the finite correlation time of the 
velocity field. 

In [[M|,[J5|], a generalized phenomenological model was considered in which the temporal 
correlation of the advecting field was set by eddy turnover (see also an earlier work [36], where 
the probability distribution function in an analogous model was studied). It was argued that 
the anomalous exponents may depend on more details of the velocity statistics, than only 
the exponents. This idea has received some analytical support in [[37]], where the case of short 
but finite correlation time was considered for the special case of a local turnover exponent. 
The anomalous exponents were calculated within the perturbation theory with respect to 
the small correlation time, with Kraichnan’s rapid-change model taken as the zeroth order 
approximation. The exponents obtained in 0 appear to be nonuniversal, through the 
dependence on the correlation time. The exact inequalities obtained in J3B] using the so- 
called refined similarity relations also point up some significant differences between the zero 
and finite correlation-time problems. 

In the paper |32|], the field theoretic renormalization group (RG) and operator product 
expansion (OPE) were applied to the model [|TUJ. The feature specific to the theory of 


turbulence is the existence in the corresponding field theoretical models of the composite 
operators with negative scaling (“critical”) dimensions. Such operators are termed “danger¬ 
ous,” because their contributions to the OPE for the structure functions and various pair 
correlators give rise the anomalous scaling, i.e., singular dependence on the IR scale with 
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nonlinear anomalous exponents. The latter are determined by the critical dimensions of 
these operators.^] The OPE and the concept of dangerous operators in the stochastic hy¬ 
drodynamics were introduced and investigated in detail in [fTf^.fTOH ; see also the review paper 

0 


and the book 42 


The part of the formal expansion parameter in the RG approach is played by the ex¬ 
ponent ( entering into the velocity correlator; see Eq. ( |1.9|) in Sec. |H| (in Ref. [|32||, it was 
denoted by e, in order to emphasize the analogy with Wilson’s e expansion). The anomalous 
exponents were calculated in |32|] to the order ( 2 of the expansion in ( for any space dimen¬ 
sion, and they are in agreement with the first-order results obtained within the zero-mode 
approach in p~^-pT| . In [|33| , the RG method was generalized to the case of a nonsolenoidal 
(“compressible”) velocity field. 

The main advantage of the RG approach (apart from its calculational efficiency) is the 
universality: it is not related to the aforementioned solvability of the rapid-change model 
and can equally be applied to the case of finite correlation time, provided the corresponding 
model possesses the RG symmetry. 


In 32 


the results were presented for the opposite 
limiting case of the time-independent (“frozen”) velocity field. 

In this paper, we apply the RG and OPE technique to the problem of a passive scalar 
held advected by a self-similar synthetic Gaussian velocity field with finite correlation time; 
the steady state is maintained by an imposed linear mean gradient. The velocity held 
satisfies a linear stochastic equation with effective viscosity and stirring force. The model 
was proposed and studied in detail (using numerical simulations, in two dimensions) in 
3; its rapid-change version is discussed in 


20.24 


We consider the problem in an 
arbitrary space dimension, d > 2; we also stress that the correlation time is not supposed 
to be small. We establish the existence in the inertial-convective range of several different 
scaling regimes and show that for some of them the structure functions and other correlation 
functions of the problem exhibit anomalous scaling behavior; we derive explicit analytical 
expressions for the corresponding anomalous exponents. 

The advection of a passive scalar held in the presence of an imposed linear gradient is 
described by the equation 


V t 0 = u 0 d 2 d - h • v, \/ t = d t + Vidi. (1.1) 

Here 6(x) = 9(t , x) is the random (fluctuation) part of the total scalar held Q(x) = 6 , (x)+h-x, 
h is a constant vector that determines distinguished direction, u 0 is the molecular diffusivity 
coefficient, d t = d/dt, <9j = d/dxi, d 2 = didi is the Laplace operator, and v(a:) = (a?)} is 

the transverse (owing to the incompressibility) velocity held. 

The velocity obeys the linear stochastic equation, cf. |J 

d t Vi + Rvi = fi , ( 1 . 2 ) 


1 For the rapid-change model, the relationship between the anomalous exponents and dimensions 

0 . 


of composite operators was anticipated in 


within certain phenomenological formulation 


of the OPE, the so-called “additive fusion rules,” typical to the models with multifractal behavior; 
see also 


4j. The RG analysis of Ref. |3^] shows that such fusion rules are indeed obeyed by 


the powers of the local dissipation rate in the model 0, and all these operators are dangerous. 
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where R [in the momentum representation R = R{k )] is a linear operation to be specified 
below and fi is an external random stirring force with zero mean and the correlator 


(. fi(x)fj(x ')) = j ^ J D f (uj, k ) exp[—i(f - t') + ik • (x - x')]. (1.3) 

Here Pij{ k) = Sij — kikj/k 2 is the transverse projector, k = |k| is the wavenumber, and 
d is the dimensionality of the x space. Following ||, we choose the correlator D * to be 
independent of the frequency, so that Eq. O contains the delta-function in time. More 
specihc, we choose 


k ) 


3 4—d—e—ri 

9o"o 


R(k ) = UqPq a k 


2 -v 


(1.4) 


where 




= y/k* 


+ m- 


(1.5) 


The positive amplitude factors g 0 (a formal small parameter of the ordinary perturbation 
theory) and w 0 are the analogs of the the coupling constant (“charge”) A 0 in the standard 
Aq 0 4 model of critical behavior, see, e.g., |[4ol.|4G|| ; in what follows we shall also term these 


parameters “coupling constants.” The exponents £ and g are the analogs of the RG expansion 
parameter e — 4 — d in the Ao0 4 model, and we shall use the traditional term “e expansion” 
in our model for the double expansion in the £-rj plane around the origin £ — r) — 0, with 
the additional convention that £ = 0{rf). The infrared (IR) regularization is provided by 
the integral scale L = 1/m; its precise form is not essential. For k » m the functions 
(1.4) take on simple powerlike form. Dimensionality considerations show that the charges 
are related to the characteristic ultraviolet (UV) momentum scale A by 


9o - A'+’j u 0 ~ A". 


( 1 . 6 ) 


From Eqs. 0) and (|T3|) it follows that v(x) obeys Gaussian distribution with zero 
mean and correlator (dropping the transverse projector) 


D v (u,k ) = 


Df(k) 


3 A—d—e—ri 

9o"o cr k 


J 1 + R 2 (k) up + [u 0 u 0 <r k r '] 2 


(1.7) 


Therefore, the exponent £ describes the inertial-range behavior of the equal-time velocity 
correlator or, equivalently, the energy spectrum 


E(k)~k d 1 J dtoD v {u, k) ~ {go^o /uq) k 1 £ , ( 1 . 8 ) 

cf. uhl where a close family of models for the velocity held has been considered for a 
strongly anisotropic shear how. The second exponent, rj, is related to the function R(k), 
the reciprocal of the correlation time at the wavenumber k (r) = 2 — z in the notation of 
|7j p|. |37l . |38|| ; our exponents are dehned so that £ = rj = 0 correspond to the starting point of 
the RG expansion). It then follows that £ = 8/3 gives the Kolmogorow “five-thirds law” for 
the spatial velocity statistics, and rj — 4/3 corresponds to the Kolmogorov frequency. 
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It was pointed out in [[3] that the linear model (|1.2|) suffers from the lack of Galilean 
invariance and therefore does not take into account the self-advection of turbulent eddies. 
It is well known that the different-time correlations of the Eulerian velocity field are not 
self-similar, as a result of these “sweeping effects,” and depend substantially on the integral 
scale; see, e.g., [|47|]. Nevertheless, the results of [[J show that the model gives reasonable 
description of the passive advection in an appropriate frame, where the mean velocity field 
vanishes. To justify the model ( |1.2|) , we also note that we shall be interested preferably in the 
equal-time, Galilean invariant quantities (structure functions, correlations of the dissipation 
rate etc.), which are not affected by the sweeping effects, and we expect that their absence 
from the model ( | 1 . 2 | ) is not essential. 

We also note that the model contains two special cases that possess some interest on 
their own. In the limit Uq —> oo, g' 0 = go/u^ = const we arrive at the rapid-change model: 


D v (u, k) -> g' 0 u 0 (k 2 + m 2 ) d/2 c/2 , C = £ 


V, 


(1.9) 


and the limit u 0 —*• 0 , = go/u 0 = const corresponds to the case of a frozen velocity field: 

D v (uj,k) —> < 7 qVo (k 2 + m 2 )~ d ^ +1 ~ £ / 2 n 8(u), (1.10) 

when the velocity correlator is independent of the time variable t — t! in the t representation. 
The latter case for h = 0 has a close formal resemblance with the well-known models of the 


random walks in random environment with long-range correlations; see [48,49 


In Sec. |TI]. we give the field theoretic formulation of the problem and discuss some its 
consequences; we also explain briefly why the ordinary perturbation theory fails to give 
correct IR behavior for some values of £ and r) and establish the relationship between the 
IR and UV problems. In Sec. E3 we discuss the UV renormalization of the model, derive 
the RG equations and present the one-loop expressions for the basic RG functions (beta 


functions and anomalous dimensions). In Sec. [TV], the analysis of the scaling behavior is 
given. Depending on the values of the exponents £ and rj entering into the velocity correlator, 
the model exhibits various types of IR scaling regimes, associated with the IR stable fixed 
points of the RG equations: 

(i) The anomalous scaling behavior with universal (in the above sense) exponents, char¬ 
acteristic of the rapid-change model, takes place for rj < e < 2 g. The anomalous exponents 
depend on the only exponent ( entering into Eq. (03). 

(ii) The anomalous scaling behavior with the universal exponents, characteristic of the 
model with time-independent (frozen) velocity field, emerges in the region e > 0 , e > 2g. 
The exponents are determined solely by the equal-time velocity correlator and depend on 
the only exponent e entering into Eq. ( |1.10|) . 

(iii) The intermediate regime with nonuniversal exponents, which depend on the am¬ 
plitudes entering into the velocity correlator, emerges for £ = 2 77 ; the Kolmogorov-type 
synthetic velocity field 0 and the case of a local turnover exponent |37[] correspond to this 
regime. The nonuniversality of the exponents in this regime is in agreement with the findings 
of Ref. [[U]], where the large d limit has been considered. [However, the exponents in our 
model turn out to be universal in the one-loop approximation], 

(iv) The diffusive-type regimes, for which the advection (i.e., the nonlinearity in Eq. 
( |1.1|) ) can be treated within the ordinary perturbation theory. These regimes take place in 
the region specified by the inequalities g > 0, g > e and g < 0, e < 0. 
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To avoid possible misunderstandings we emphasize that the limits g 0j u 0 — 1 ► 0 or g Q , 
u 0 —■► oo are not supposed to be performed in the original correlation function (|1.7|) ; the 
parameters g 0j u 0 are fixed at some finite values. The behavior specific to the models (p..9|), 
( |1.10| ) arises asymptotically in the regimes (i) and (ii) as a result of the solution of the RG 
equations, when the “RG flow” approaches the corresponding fixed point. Therefore, we 
deal with the finite correlation time, and there is no problem with the steady state in the 
frozen case even in two dimensions. The regions of IR stability of the regimes (i)-(iv) in the 
£-rj plane, given above, are identified to the first order of the e expansion, but some of their 
boundaries are found exactly. 

In the regimes (i)-(iii), the correlation functions of the model exhibit anomalous scaling 
behavior, i.e., singular dependence on the IR scale m with nonlinear “anomalous exponents.” 
Within the RG and OPE approach, the latter are related to the scaling dimensions of the 
tensor composite operators dO ■ ■ ■ d9\ these dimensions are calculated explicitly to the first 
order of the e expansion (one-loop approximation) in Sec. [V]. The inertial-convective-range 
asymptotic expressions for the structure functions of arbitrary order (even and odd) and the 
equal-time correlations of the scalar gradients are obtained in Sec. Q using the OPE. 

As the exponents e and r\ increase, the powers of the velocity held also become danger¬ 
ous, and their contributions to the OPE should be summed. The required summation is 
performed in Sec. ED on the example of the second-order structure function in the “frozen” 
regime; for the rapid-change regime the problem is absent. This summation might be inter¬ 
esting as a possible model of the origin of the anomalous scaling in the structure functions 
of the velocity itself: it was argued in J5(| that the singular m dependence of the equal-time 
correlators for the stochastic Navier-Stokes (NS) equation is related to infinite families of 
dangerous operators. 

The formulation ( |1.6| ) is typical to the models of critical behavior [i.e., the dimensional 
coupling constants are expressed only through the UV scale], and we shall call it “standard.” 
In d|. |?8| . pd] |o8l, a different version of the problem was considered, in which the velocity 
correlator contains nontrivial dependence on the integral turbulence scale. This “exotic” 
(from the viewpoints of the theory of critical behavior) formulation requires special attention; 
it is considered briefly in Sec. 01 

The results obtained are reviewed in Sec. 0 where we also discuss briefly the pas¬ 
sive advection by the non-Gaussian velocity field governed by the nonlinear stochastic NS 
equation. Our approach is generalized directly to this case, and the explicit expressions for 
the anomalous exponents can readily be obtained in the first order of the corresponding e 
expansion. We also discuss new problems that arise in the NS model beyond the £ expansion. 


II. FIELD THEORETIC FORMULATION OF THE MODEL. IR AND UV 
SINGULARITIES IN THE PERTURBATION THEORY 


According to the general theorem (see, e.g., Refs. [[f5j [f6|), the stochastic problem (|1.1|) ~ 
( |1.3| ) is equivalent to the field theoretic model of the doubled set of fields = {9, 9f v, v'} 
with action functional 


S($) = (1/2 )v'DV + v'[-<9 t v - Rv] + 9' \-d t 9 - {wd)9 + v 0 d 2 6 - h • v 


( 2 . 1 ) 
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Here Df is the correlator 0, the required integrations over x = (t, x) and summations 
over the vector indices in Eq. (2.1) and analogous formulas below are implied. 

The formulation (|2.1| ) means that statistical averages of random quantities in the stochas¬ 
tic problem (|1.1|) - (|1.3|) coincide with functional averages with the weight exp S'(<£>), so that 
generating functionals of total [G(H)] and connected [bE(v4)] Green functions are represented 
by the functional integral 


G(A) = exp W(A) = j P$exp[S($) + A$\ 
with arbitrary sources A(x) in the linear form 

H<f> = f dx[A e (x)6(x) + A d \x)6'(x) + AJ(x)vi(x) + A^' (x)v((;x)\. 


( 2 , 2 ) 


(2,3) 


In the following, we shall not be interested in the Green functions involving the auxiliary 
vector field v', so that we can set A v = 0 in Eq. ( ]2.3| ) . It is then convenient to perform the 
Gaussian integration over v' in Eq. (|2.2| ) explicitly. We arrive at the field theoretic model 
of the reduced set of fields = {9, 6 ', v} with the action 


S($) = d' 


—d t 6 — (vd)9 + v 0 d 2 9 — h • v 


vD„-‘v/2. 


(2,4) 


The first four terms in Eq. (|2.4|) represent the Martin-Siggia-Rose-type action for the 
stochastic problem (TTT ) at fixed v, and the last term represents the Gaussian averaging 
over v with the correlator D v from Eq. ( |1 . 7| ) . 

The model ( |2.4| ) corresponds to a standard Feynman diagrammatic technique with the 
triple vertex —6'(yd)6 = 6'VjVjO with vertex factor 


Vj = i kj, (2.5) 

where k is the momentum flowing into the vertex via the field 9', and the bare propagators 
(in the momentum-frequency representation) 

(66') 0 = (6'ey 0 = (-iu + v 0 k 2 )- 1 , 

(99) o = (99') 0 h i h j (v i v j )o(9 , 9)o, 

(Ovi) o = -(99') 0 h j (v j v i ) 0 , 

(9’9') 0 = 0, (2.6) 

where hi is a component of the vector h and the bare propagator (ihVj) o is given by Eq. 

0 , 

The magnitude h = |h| can be eliminated from the action (|2.4|) by rescaling of the 
scalar fields: 9 —> h9 : 9' —> 9'/h. Therefore, any total or connected Green function of 
the form (9(xi) ■ ■ ■ 9(x n ) 9'(y\) ■ ■ ■ 9'(y p )) contains the factor of h n ~ p . The parameter h 
appears in the bare propagators ( |2.6| ) only in the numerators. It then follows that the 
Green functions with n — p < 0 vanish identically. On the contrary, the 1-irreducible func¬ 
tion (9(x i) • ■ • 9(x n ) 9'(yi) ■ ■ ■ 9'(y p ))i_ h . contains a factor of h p ~ n and therefore vanishes for 
n — p > 0 ; this fact will be relevant in the analysis of the renormalizability of the model (see 
below). 
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Another important consequence of the representation (|2.2|) , (|2.4| ) is that the large-scale 
anisotropy persists, through the dependence on h, for all ranges of momenta (including con¬ 
vective and dissipative ranges), and that the dimensionless ratios of the structure functions 
are strictly independent on h ; cf. p[-|5|,|34|-06l] . It is noteworthy that all these statements 
equally hold for any statistics of the velocity held (not necessarily Gaussian or synthetic), 
provided its distribution is independent of h. 

However, the ordinary perturbation theory fails to give correct IR behavior of Green 
functions for some values of the exponents £ and 77 . This can easily be illustrated on the 
simplest example of the 1-irreducible Green function {O'9) i_ ir . It satisfies the Dyson equation 
of the form 


{O'O) i_ ir = -iu + u 0 k 2 - E 0 / 0 (a;, k), (2.7) 

where Eg'e is the self-energy operator represented by the corresponding 1 -irreducible dia¬ 
grams. Its one-loop approximation has the form 

. ( 2 - 8 ) 


Here and below the solid lines in the diagrams denote the bare propagator {66') o from Eq. 
( pT 6 [) , the end with a slash corresponds to the held 6', and the end without a slash corresponds 
to 9 ; the dashed lines denote the bare propagator (|E7|); the vertices correspond to the factor 
( j2.5|) . The analytic expression for the diagram in (|2.8|) has the form 


Ve(o;,A;) = -kjkj J J 


dq 


Fiji q) D v {uj',q ) 


{2n) d -i(u; + u/) + i/ 0 (q + k ) 2 ’ 


(2.9) 


where q = |q| and D(iv',q) is given by Eq. (|1~7|); the factor of kikj arises from the vertex 
factors fl2.5|) . Integration over u/ in Eq. (|2.9| ) yields 


S 0 ' 6 »(u;, k) 


k k 90*0 f dq _ Pjj{q)cT 2 q d £ 

' j 2 u 0 J {2-kY —\tjj + z/ 0 (q + k) 2 + u 0 u 0 al~ V 


( 2 . 10 ) 


We are interested in the IR behavior of the function ( |2.10| ), i.e., the behavior of small k, 
iv and m. It is easily seen that this behavior is nontrivial in the region on the e-r) plane, 
determined by the inequalities 77 < 0 , e > 0 and 77 > 0 , e > 77, because the integral in (|2.10|) 
is then IR divergent if k, iv and m are simply set equal to zero. On the contrary, for the 
rest of the e-r) plane, the leading term of the desired asymptotic behavior is indeed obtained 
simply by setting k = uj — m — 0. The analysis is extended directly to the higher-order 
diagrams; it shows that these IR singularities enhance as the order of a diagram increases, 
and that they take place only within the same region on the e-rj plane. The IR singularities 
compensate the smallness of the coupling constant go, assumed within the framework of the 
ordinary perturbation theory. Therefore, in order to find correct IR behavior we would have 
to sum the entire series even if the expansion parameter, go, were small. 

It is also clear that these IR singularities get weaker as the parameters e, 77 decrease, and 
they would disappear at e = g = 0 if we could take this limit in Eq. ( |2.10| ). However, this 
is impossible owing to the UV divergence in the integral ( |2.10|) at this point. I 11 general, 
the diagrams of are UV divergent in the region 77 > 0, e < 0 and 77 < 0, e < 77, and the 
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UV cutoff at q = |q| ~ A is then implied in (|2.10|) and higher-order diagrams. If the point 
e — rj — 0 is approached from inside the region of UV convergence, the UV singularities 
manifest themselves as poles in e, ij and their linear combinations. The elimination of these 
poles is the classical UV problem, and its solution is given by the standard theory of UV 
renormalization; the RG equations are obtained within the framework of this theory and 
express the simple idea of nonuniqueness of the renormalization procedure. The correlation 
between the IR and UV singularities near the “logarithmic point” e = rj = 0, noted above, 
explains to some extent why the RG method, which is closely related to the UV divergences, 
can be a useful tool in studying the IR behavior, and why the exponents £ and r] are expected 
to be relevant small parameters in the RG expansions. 

Surprisingly, simple arguments given above lead to reasonable conclusions: the rigorous 
RG analysis confirms that the Green functions of the model indeed show anomalous IR 
behavior for some values of £ and 77 , and the region determined by the inequalities 7] < 0 , 
£ > 0 and r) > 0, e > rj coincides with the region of stability of the corresponding fixed 
points in the linear approximation; see Sec. 


III. IV. 


III. UV RENORMALIZATION OF THE MODEL. RG FUNCTIONS AND RG 

EQUATIONS 


The renormalization of the model (|2.4j) is similar to the renormalization of the simpler 
rapid-change model, considered in detail in |32[; below we confine ourselves to only the 
necessary information. 

The analysis of UV divergences is based on the analysis of canonical dimensions, see 
46| , [5l|| . Dynamical models of the type (|2.4| ), in contrast to static models, are two-scale 
4l| , |42|j5^| , i.e., the action functional 


is invariant with respect to the two independent 
scale transformations, S(Q',z'i) = S($,Zi), where $ = {9,9 ', v} and Z{ = {g 0 , uq, v 0 , m} is 
the full set of the model parameters. In the first transformation, the time variable is fixed 
and the space variable is dilated along with all the fields and parameters: 


x) —> <h'(t,x) = AAx), Zi —> z[ = 


dt 


(3.1) 


and in the second the space variable is fixed and all the other quantities are dilated: 

<h(t,x) —► &(t, x) = A d *<f>(Af,x), Zi —> z'i = A^*» Zi . (3.2) 


Here A > 0 is an arbitrary transformation parameter, and two independent canonical dimen¬ 
sions, the momentum dimension dp and the frequency dimension dip, are assigned to each 
quantity F (a field or a parameter in the action functional). These canonical (“engineering”) 
dimensions should not be confused with the exact critical dimensions: the latter are subject 
to nontrivial calculation, while the former are simply determined from the natural normal¬ 
ization conditions d\ = — — 1 , d% — d“ — 0 , d 1 ^ = d{ = 0 , d “ = —= 1 , and from the 
requirement that each term of the action functional be dimensionless [i.e., be invariant with 
respect to the transformations (0) and (|3.2|) separately]. Then, based on dp and dp, one 
can introduce the total canonical dimension pT|JT 2 ||) 2 fl , which corresponds to the dilatation 
with fixed value of u 0 (i.e., zero canonical dimension can be assigned to u 0 ). In our model, 
d t oc u 0 d 2 , so that the total dimension is given by dp = d k F + 2 dp. 
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In the action (|2.4|) , there are fewer terms than fields and parameters, and the canonical 
dimensions are not determined unambiguously. This is of course a manifestation of the fact 
that the “superfluous” parameter h — |h| can be eliminated from the action; see above. 
After it has been eliminated (or, equivalently, zero canonical dimensions have been assigned 
to it), the definite canonical dimensions can be assigned to the other quantities. They are 
given in Table |. including the dimensions of renormalized parameters, which will appear 
later on. 

From Table | it follows that the model is logarithmic (the both coupling constants go 
and Uo are dimensionless) at e = rj = 0. This means that the UV divergences in the Green 
functions have the form of the poles in e, p, and all their possible linear combinations. 

The total dimension dp plays in the theory of renormalization of dynamical models the 
same role as does the conventional (momentum) dimension in static problems. The canonical 
dimensions of an arbitrary 1-irreducible Green function T = (<!>... $)!_ ir are given by the 
relations 


4 = d - N*d%, 


d£ = 1 - N^d%, 


dp — d r + 2 dp = d + 2 — N$d, 


■J)U(j>, 


(3,3) 


where N$ = { Ng , Ng>, N v } are the numbers of corresponding fields entering into the function 
T, and the summation over all types of the fields is implied. The total dimension dr is 
the formal index of the UV divergence. Superficial UV divergences, whose removal requires 
counterterms, can be present only in those functions T for which dp is a non-negative integer. 

Analysis of divergences in the problem (|2.4| ) should be based on the following auxiliary 
considerations; cf. 0.0.0: 

(1) All the 1-irreducible Green functions with Ngi < Ng vanish; see Sec. pQj. 

(2) If for some reason a number of external momenta occur as an overall factor in all the 
diagrams of a given Green function, the real index of divergence d' v is smaller than dr by the 
corresponding number (the Green function requires counterterms only if d' r is a non-negative 
integer). 

In the model (|2.4| ), the derivative 8 at the vertex Q'(yd)6 can be moved onto the held 
O' by virtue of the transversality of the held v. Therefore, in any 1-irreducible diagram 
it is always possible to move the derivative onto any of the external “tails” 9 or 9', which 
decreases the real index of divergence: d' T = dr — Ng — Ng>. This also means that the helds 
9, 9' enter into the counterterms only in the form of the derivatives 89 and 89'. 

From the dimensions in Table | we find dp = d + 2 — N v + Ng — (d + l)Ngi and d' r = 
(d+2)(l — Ngr) — N v . From these expressions it follows that for any d, superhcial divergences 
can exist only in the 1-irreducible functions (O'6... 9) i_; r with Ngr = 1 and arbitrary value 
of Ng , for which dp = 1 + Ng, d' T = 0. However, all the functions with Ng > Ngr vanish (see 
above) and obviously do not require counterterms. As in the case of the rapid-change model 
|32|,|33|1, we are left with the only superficially divergent function (9'9)i_ n \ the corresponding 
counterterm contains two symbols 8 and is therefore reduced to 9'8 2 9. The inclusion of this 
counterterm is reproduced by the multiplicative renormalization of the parameters go, uq, 
and Vq in the action functional (|2.4|): 


Vo = vZ v , 


= qu £+v Z, 


9 o = 99 


gi 


u 0 ug' Z u , 


(3.4) 


where the dimensionless parameters g , u, and v are the renormalized analogs of the bare 
parameters, /i is the renormalization mass in the minimal subtraction (MS) scheme, which 
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we always use in practical calculations, and Zi = Zi(g,u ) are the renormalization constants. 
They satisfy the identities 


Z 9 = ^ 7 3 , 


Z u — z. 


-i 


(3.5) 


which result from the absence of the renormalization of the contribution with D„ in the 


functional (2.4). No renormalization of the fields, the “mass” m, and the vector h is required, 
i.e., = 1 for all <f> and Z m — Z h — 1. 

The renormalized action functional has the form 


S' ren ($) = 9' -d t 9 - (’ vd)9 + vZ v d 2 9 - h v - vD~W/ 2, 


(3.6) 


where the correlator D v is expressed in renormalized parameters using the formulas fl3.4Q : 


D v (u, k ) 


3 4 —d—e—rt 

gv 6 g, +n a k 
u 2 + [uv n n crl^ v ] 2 


(3.7) 


The relation S(&,eo) = S' ren ( < h, e, g) (where e 0 is the complete set of bare parameters, 
and e is the set of renormalized parameters) for the generating functional W (A) in Eq. (|2.2|) 
yields W(A, e 0 ) = W ren (A, e, /i). We use V tl to denote the differential operation for hxed 
eo and operate on both sides of this equation with it. This gives the basic RG differential 
equation: 


TArg W ien (A, e, n) — 0, 

where T>rg is the operation T> fI expressed in the renormalized variables: 


(3.8) 


Vrg = V,,. + p g (g, u)d g + {3 u (g, u)d u - 7 u (g, u)T> v . (3.9) 

In Eq. ( |3.9D , we have written D x = xd x for any variable x, and the RG functions (the (3 
functions and the anomalous dimension 7 ) are defined as 

7 „ = X> At ln Z„, (3.10a) 

^ = ^ = ^[-£-^ + 37,,], (3.10b) 

(3 U = V^u = 4-77 + 7 „]. (3.10c) 


The relations between (3 and 7 in Eq. ( 3.10 ) result from the definitions and the relation 

(PD- 

Now let us turn to the explicit calculation of the constant Z v in the one-loop approx¬ 
imation in the MS scheme. The constant Z v is determined by the requirement that the 
1-irreducible Green function (9'9) i_ ir , when expressed in renormalized variables, be UV fi¬ 
nite [i.e., have no singularities for e, 77 —*■ 0]. The Dyson equation (2.7) relates this function 
to the self-energy operator H$ig, and Eq. ( |2.10| ) gives the explicit expression for the latter in 
the first order O(go) of the unrenormalized perturbation theory. Now we have to calculate 
the function in the order O(g) of the renormalized perturbation theory; therefore we 
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should simply replace u 0 —> ► v in the propagator (69') 0 and use the expression ( |3.7| ) for the 
velocity correlator in Eq. (|2.8|) , which leads to the substitution g 0 —> g/J. £+v , u 0 —> w/T 7 , 
Uq — > v in Eqs. ( |2.9| ), (|2.10| ). We know that the divergent part of the diagram is indepen¬ 
dent of u, so that we can set oj = 0 in what follows. It is also convenient to cut off the 
integral over q from below at q ~ m and set m — 0 in the integrand (the integral diverges 
logarithmically, and its UV divergent part is independent of the specific form of the IR regu¬ 
larization). Furthermore, we can set k = 0 in the integrand (we know that the counterterm 
is proportional to k 2 , and the factor of k 2 has already been isolated from the integral) and 
make use of the isotropy, namely, 

f dqf(q)Pij(q ) = j /(<?)• 


Then Eq. Q2.10|) yields 


E e / 0 (o; = 0, k) ~ —k 


2 fP J l j£ (d ~ 1 ) J 


2nd 


(3.11) 


where we have written 


J = 


dq q d £ 

( 2n) d 1 + u (/d/q)^ 


(3.12) 


and ~ denotes the equality up to the UV finite parts. The expansion of the integrand in u 
gives 


J = £(-“)> 


ST] 


*> 


s=0 


(2t) 


e—sr) 


n oo 


, i_i sri m £ sri 

£ + ST] 


(3.13) 


where the parameter m arises from the IR limit in the integral over q and Sd = 2ir d P/T(d,/2) 
is the surface area of the unit sphere in d -dimensional space. 

Finally, from Eqs. (|3.11|) and (|3.13|) we obtain 


„ , _ , N - agvk 2 

E e / 0 (a; = 0, k) ~- J2 ' 

u s =o 


-u) s (g,/m) £+sri 

£ + ST] 


(3.14) 


where we have written 


(d-l)S d 

2d(2n) d ‘ 


(3.15) 


The renormalization constant Z v is found from the requirement that the UV divergences 
cancel out in Eq. (EZP after the substitution u 0 = vZ v . This determines Z v up to an UV 
finite contribution; the latter is fixed by the choice of the renormalization scheme. In the 
MS scheme all the renormalization constants have the form “1 + only poles in £, 7] and their 
linear combinations,” which gives the following expression 


ag_y, (~vy_ 

u “ £ + si] ’ 


(3.16) 
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with the coefficient a from Eq. (|3.15|) . 

In contrast to the rapid-change model, the one-loop approximation in the case at hand 
is not exact: the expression (|3.16|) has nontrivial corrections of order g 2 , g 3 , and so on. 
The series in Eq. ( |3.16|) can be expressed in the form of a single integral, but this is not 


convenient for the calculation of the RG functions. 

The RG functions in the one-loop approximation can be calculated from the renor¬ 
malization constant ( f3.16| ) using the identity V M = j3 g d g + /3 u d u , which follows from the 
definitions (|3.10|) and the fact that Z u depends only on the charges g, u. Within our accuracy 
this identity is reduced to ~ — (e + g)Z~> g — gT> u . From Eq. (fj.ltjp it then follows: 


lv = 


{e + rj)V g + rjV u — ]T 


QjQ ^( ZL ) 

u + 




ag 


(3.17) 


up to the corrections of order g 2 and higher. 
( P.l?l ) using the relations ( |3.10b|) , (|3.10c|) - 


The beta functions are obtained from Eq. 


IV. FIXED POINTS AND SCALING REGIMES 


It is well known that possible scaling regimes of a renormalizable model are associated 

HI 


The 


with the IR stable fixed points of the corresponding RG equations, see, e.g., 
fixed points are determined from the requirement that all the beta functions of the model 
vanish. In our model the coordinates g*,w* of the fixed points are found from the equations 


/ %{g *, «*) = Pu{g*, u*) = 0 


(4.1) 


with the beta functions given in Eqs. (|3.10b| ), (|3.10c|) . The type of the fixed point is 
determined by the eigenvalues of the matrix = dfli/dgi}, where Pi denotes the full 

set of the beta functions and g, is the full set of charges. For the standard (as in Eq. ([Q|)) 
formulation of the problem the IR asymptotic behavior is governed by the IR stable fixed 
points, i.e., those for which all the eigenvalues are positive. 

From the equations ( |3.10b|) , (|3. 10c|) we obtain the exact relation j3 g / g — 3(3 u /u = 2g — e. 
It shows that the beta functions (3 g , /3 U cannot vanish simultaneously for finite values of 
their arguments. [The only exception is the case 2 r) = e. We shall study it separately, and 
for now we assume 2r/ ^ £.] Therefore, to fold the fixed points we must set either u = 0 or 
u = oo and simultaneously rescale g so that the anomalous dimension remain finite. 

In order to study the limit u —> oo we change to the new variables w = 1/u, g' = g/u 2 ; 
the corresponding beta functions have the form 


Pw = V g w = -f3 u /u 2 = w[r)~ 7 „], 
Pg' = V g g' = pg/u 2 - 2 gPujuP = g'[rj - £ + 7 „], 


and for the one-loop anomalous dimension we obtain from Eq. ( |3.17|) 

lv = ag'/{l + w) 


(4.2) 


(4.3) 


with the constant a defined in Eq. ( |3.15|) . From the expressions (|4.2| ) we hnd two fixed 
points, which we denote FPI and FPII. The first point is trivial, 
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FPI : 


w* = g* = 0 ; 


7 ^ = 0 . 


(4.4) 


The corresponding matrix is diagonal with the diagonal elements 

Oi = r /, 0.2 = f] ~ £■ (4.5) 

For the second point we obtain 

FPII : w* = 0, g* = (e-rj)/a\ 7 t = £ ~ V- (4.6) 

The corresponding matrix hi is triangular, d g >fJ w = 0, and its eigenvalues coincide with the 
diagonal elements: 


= d w p w = 77 - 7 * = 2r] - e, 
0 2 = dg'fig' = ag'* = e - 77 . 


(4.7) 


We note that the expressions for 7 * in Eq. ( |4.6| ) and for ffi in ( |4.7| ) are exact, i.e., they have 
no corrections of order 0 (e 2 ) [we take e ~ 7 , so that here and below 0(s 2 ) denotes all the 
terms of the form eg, rj 2 and higher]. 

Now let us turn to the regime with u —> 0. In order to study this limit we change to the 
new variable g" = gju\ the corresponding beta functions have the form 


Pg» = V^g" = Pg/u - gp u /u 2 = g"[-e + 2^ v ], 

Pu = u[—rj + 7 ^] 


(4.8) 


[the function P u is the same as in Eq. ( p,10c|) ]. The one-loop anomalous dimension ( |3.17|) 
takes the form 


lu = ag"/(l + u). 


(4.9) 


From the expressions (T 8 ) we find two fixed points, which we denote FPIII and FPIV. The 
first point is trivial, 

FPIII: m* — g” — 0; 7 ^ = 0. (4.10) 

The corresponding matrix 0 is diagonal with the elements 

Ill = —e, 0 2 = —TJ. (4-11) 

For the nontrivial point we obtain 

FPIV: u, = 0, g'l = e/2a; H = e/2. (4.12) 

The corresponding matrix 0 is triangular, d g nP u = 0, and its eigenvalues have the form 


fii = duPu = -7 + 7 * = (e - 2rj)/2, 
0 2 = dg"P g " = 2 ag” = e. 


(4.13) 


The expressions for 7 * in Eq. (|4 .1 2|) and for in Eq. ([4.13|) are exact. Of course, the 
expressions ( |4.5|) , ( |4.11|) , and 7 * = 0 for the trivial fixed points are also exact. 
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In the special case e = 2r/ the beta functions (|3.10b| ), (|3.10c|) become proportional, and 
the set (Hi reduces to a single equation. As a result, the corresponding nontrivial fixed 
point, which we denote FPV, is degenerate: rather than a point, we have a line of fixed 
points in the g-u plane. It is given by the relation 

FPV : g*/u*(u* + 1) = 77 /a; 7 * = g = e/2. (4.14) 

The exact expression for 7 * follows from the relation between the RG functions in Eq. ( [3.10D . 
The eigenvalues of the matrix (which is not diagonal here) have the form 

= 0, 0,2 — V (2 + u*)/(l + u*). (4.15) 


The vanishing of the element IR reflects the existence of a marginal direction in the g-u plane 
(along the line of the fixed points) and is therefore an exact fact. The coordinates of a point 


on the line (|4.14|) can also be expressed explicitly as functions of the dimensionless parameter 
p = go/uQ using the exact relation go/u^ = g*/u\. The actual expansion parameter appears 
to be y/rj rather than 1 ] itself, and the zeroth order approximation has the form 


g* = ( v/a) 3/2 p 1/2 , w* = (rj/ ap) 1/2 , 0 2 = 2 rp 


(4.16) 


I 11 Figure I, we show the regions of stability for the fixed points FPI-FPV in the e-ij 
plane, i.e., the regions for which the eigenvalues of the matrix are positive. The boundaries 
of the regions are depicted by thick lines. We note that the regions adjoin each other without 
overlaps or gaps. This fact is exact for the ray e = 2rj > 0, the boundary between the regions 
of stability for the points FPII and FPIV [at the same time, this ray is the region of stability 
for the point FPV], On the contrary, the boundary e = 77 , r\ > 0 for the point FPII and £ — 0 , 
e > 7] for FPIV are approximate, so that the gaps or overlaps can appear in the two-loop 
approximation. The regions denoted as FPIVa and FPIV 6 with the boundary £ = 2 both 
correspond to the same fixed point FPIV; the part FPIV 6 represents the region in which 
the velocity held has negative critical dimension; see Sec. E3 

Surprisingly, Fig. I has some resemblance with the phase diagrams presented in Refs. 0,0, 
despite the essential difference between the models (in those papers, a strongly anisotropic 
velocity held has been studied). Indeed, the boundaries between the diffusive-type behavior 
(“homogenization regime” in terminology of J7[] ) and convective-type regimes (“superdiffu- 
sive behavior”) in the two models coincide (however, in our case they are not exact and will 
be affected by the 0(e) corrections). Furthermore, the Kolmogorov point ( e = 8/3, rj = 4/3) 
in our case and in |7j lies on a boundary between two nontrivial regimes. We also note that 
the boundary 2rj = £ between the rapid-change and frozen regimes was anticipated on phe¬ 
nomenological grounds in Ref. j37j, see also |I 8 |; their arguments can be linked directly to 
the RG analysis (see below). 

It is clear from the definition of the parameters g\ g" that the critical regime governed 
by the point FPII corresponds to the rapid-change limit (|1.9| ) of our model, while the point 
FPIV corresponds to the limit of the frozen velocity held; see Eq. (|1.10|) . This shows that in 
the latter case, the temporal fluctuations of the velocity held are asymptotically irrelevant in 
determining the inertial-range behavior of the scalar, which is then completely determined 
by the equal-time velocity statistics. In the former case, spatial and temporal fluctuations 
are both relevant, but the effective correlation time of the scalar held becomes so large under 
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renormalization that the correlation time of the velocity can be completely neglected. The 
inertial-range behavior of the scalar is determined solely by the u = 0 mode of the velocity 
held; this is the case of the rapid-change model. 

We then expect that all the critical dimensions at the point FP11 [FPIV] depend on the 
only exponent ( = £ — r] [e] that survives in the limit in question, and coincide with the 
corresponding dimensions obtained directly for the models (|1.9|) [( |1 .1 Q|) 1 . This is indeed the 
case; see Eqs. ( ff.20p and ( |5.14|) below. 

In the regimes governed by the trivial fixed points FPI and FPI1I, the contribution of 
the convection dies out in the IR asymptotic region; the 1R behavior has purely diffusive 
character, while the convection can be treated within ordinary perturbation theory. The 
existence of the two fixed points, the frozen and the rapid-change ones, implies that for 
rj < 0 transport by small wavenumbers k —» 0 is governed by equal-time (spatial) velocity 
statistics, while for 7] > 0 transport by small wavenumbers is determined by the u = 0 mode, 
i.e., the time decorrelated component of the velocity field. If there were IR singularities in 
the scalar correlations, they would be determined by the contributions of small momenta, 
and these two regimes would be really different. However, in the regions of stability of the 
trivial fixed points there are no such singularities (see the discussion in Sec. Q). Moreover, 
in these regimes all momenta k contribute to the long-term, large-scale transport properties 
of the scalar held (we recall that for rj > 0, £ < 0 and r) < 0, e < rj, the actual UV cutoff A 
has to be introduced, see Sec. p. and the main contribution to the perturbative diagrams 
then comes from the momenta of order k ~ A). The RG is not suitable for studying such 
“A divergent,” analytic in momenta and frequecies, quantities. Therefore, the splitting of 
the homogenization regime into the rapid-change and frozen parts is not meaningless, but 
not practically useful. Probably for this reason it was not mentioned in Refs. M- In what 
follows, we shall focus our attention on the nontrivial (anomalous) regimes. 

The solution of the RG equations in conformity with the stochastic hydrodynamics is 
HT)| [T3| in detail; see also [p5|.(T]|| for the case of the rapid-change models. 


discussed in Refs. 

Below we restrict ourselves with the only information we need. 

Any solution of the RG equation (|3.8| ) can be represented in terms of invariant variables 
g(k), u(k ), and P(k), i.e., the first integrals normalized at k = /i to g, u, and u, respectively 
(we recall that /i is the renormalization mass in the MS scheme). The relation between the 
bare and invariant charges has the form 


go = k e+ " g Z,l g u), u 0 = k n u Z u (g, u), = vZ„(g, u), 


(4.17) 


see, e.g., PT|J42| , |5Uf . Equation ( [4.17| ) determines implicitly the invariant variables as functions 
of the bare parameters; it is valid because both sides of it satisfy the RG equation, and 
because Eq. ( |4. 17|) at k — /j coincides with (|3.4|) owing to the normalization of the invariant 
variables. 

Correlation time of the velocity field at the wavenumber k is determined by the rela¬ 


tion t v 1 (k) = R{k ) = uou 0 k 2 see Eqs. (|1.4|), (|1.7| ). Correlation time of the free scalar 


held is given by fg 1 (/c) = u 0 k 2 , in the presence of advection it is replaced by the exact 
expression tg l (k ) = u(k)k 2 . The relations (175) and ( 4. 17|) allow the bare parameters and 
renormalization constants to be eliminated from the ratio tg(k)/t v (k ); this gives 


t e (k)/t v (k) = u{k ) oc const k ??+7l/ 


(4.18) 
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The last relation in Eq. (|4.18|) holds for k —» 0. It follows from the RG equation T> k u = 
P u (g,u), which reduces to V k u = u[—g + 7*] near a fixed point; see Eq. (|3.10c|) . Equation 
( |4.18| ) discloses the precise physical meaning of the invariant variable u: the ratio of the 
velocity and scalar correlation times at the wavenumber k. Now we can complete the above 
discussion of the scaling regimes and relate it to the phenomenological arguments given in 
Refs. ||37l , [38|| . From (|4.18|) it follows that for the fixed points FPI and FPII the velocity 
correlation time t v {k) becomes very small in comparison to tg(k) for k —> 0 and can be 
disregarded; we arrive at the time-decorrelated velocity field. For FPIII and FPIV, the 
opposite inequality, t v (k) » tg(k ), holds for small momenta, the temporal fluctuations of 
the velocity are “frozen in,” and its correlation time can be replaced with t v (k) = 00 . [Using 
the representation ( jl.l§|) and the exact expressions for 7 * in Eqs. (|4~4|), (|4.6|) , ( |4.1C| ) and 
( |4.12| ) , one can easily check that u — > 00 for FPI and FPII and u —> 0 for FPIII and FPIV, 
in agreement with the analysis of the U matrix.] However, these strong inequalities for the 
correlation times hold only asymptotically for k —> 0 , and therefore the exact correlator (O) 
can be replaced with its limits ( |1.9| ) or ( |1.1C| ) only in calculation of a quantity dominated 
by small k modes of the velocity field. Finally, for the point FPV one has 7 * = rj and the 
ratio (|4.18|) remains finite for k —> 0; this is the case of the local turnover exponent, studied 

B7I. 


m 


Let F be some multiplicatively renormalized quantity (a parameter, a held or composite 
operator), i.e., F = ZpF ren with certain renormalization constant Z F . Then its critical 
dimension is given by the expression 


A[F] = A F = dp + A^dp + 7 p, 


(4.19) 


see, e.g., 


40-42,52 . 


Here 


and dp are the corresponding canonical dimensions, 7 ^ is 


the value of the anomalous dimension 7 p(g) = T> tI In Zp at the fixed point in question, and 
A^ = 2 — 7* is the critical dimension of frequency. For the nontrivial fixed points we obtain 


An, = 2 



for FPII, 
for FPIV, 
for FPV 


(4.20) 


(we recall that ( = e — r), see (|1.9|) ). The critical dimensions of the fields $ in our model are 
also found exactly: 


A v = 1 — ll 


Aa = -1 


Agi — d + 1, 


(4.21) 


and for the IR scale we have A m = 1 [we recall that all these quantities in the model 
o are not renormalized and therefore their anomalous dimensions vanish identically, 
7 $ m = 0]. It is also not too difficult to show that the composite operator 6 n in the model 
( p.4|) is not renormalized and therefore its critical dimension is given simply by the relation 
A[d n ] = nA{6}] cf. [[42| for the rapid-change case. 

We note that the canonical dimensions of the fields 9, 6' in our model (see Table I) differ 
from their counterparts in the isotropic rapid-change model (see Table I in Ref. |J2]]). As a 
result, the critical dimensions (|4.20|) and ( 4.21|) at the point FPIV differ from their analogs 
for the rapid-change model, in spite of the fact that the anomalous dimensions are identical. 
In principle, the canonical dimensions in two models can be made equal by an appropriate 
rescaling of the scalar fields; we shall not dwell on this point here. 
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Let G[r) = (Fi(x)F 2 (a;')) be an equal-time two-point quantity, for example, the pair 
correlation function of the primary fields $ or some multiplicatively renormalizable com¬ 
posite operators. The existence of a nontrivial 1R stable fixed point implies that in the IR 
asymptotic region Ar >> 1 and any fixed mr the function G(r ) takes on the form 

G{r) ~ Vq g A dG (Ar) _A,G ^(mr), (4.22) 

with the values of the critical dimensions that correspond to the fixed point in question and 
certain scaling function £ whose explicit form is not determined by the RG equation itself. 
The canonical dimensions cIq, da and the critical dimension A^ ot the function G(r) are 
equal to the sums of the corresponding dimensions of the quantities iq. 


V. CRITICAL DIMENSIONS OF THE COMPOSITE OPERATORS 89 ■ ■ • 89 


In the following, an important role will be played by the composite operators of the form 

F[n,p] = d h 9 ■ ■ ■ d ip 9 (diOdiO) 1 , (5.1) 

where p is the number of the free vector indices and n = p + 21 is the total number of the 
fields 9 entering into the operator; the vector indices of the symbol F[n,p] are omitted. 

Coincidence of the held arguments in Green functions containing a composite operator 
F gives rise to additional UV divergences. They are removed by a special renormalization 
procedure, described in detail, e.g., in [|45|,f|(||)I|] . The discussion of the renormalization of 

mm 


composite operators in turbulence models can be found in 


see also Ref. 


for 


the case of Kraichnan’s model. Owing to the renormalization, the critical dimension A [F] 
associated with certain operator F is not in general equal to the simple sum of critical 
dimensions of the fields and derivatives entering into F. As a rule, the renormalization of 
composite operators involves mixing, i.e., an UV finite renormalized operator is a linear 
combination of unrenormalized operators, and vice versa. 

The analysis of UV divergences is related to the analysis of the corresponding canonical 
dimensions, cf. Sec. m It shows that the operators F[n,p\ mix only with each other in 
renormalization, with the multiplicative matrix renormalization of the form (dropping the 
vector indices everywhere) 


U[?T.,p] ^[n,p] [n',p']Fren \p> iP\ (5-2) 

Here F ren is the renormalized analog of the operator F and Z is the matrix of renormalization 
constants. The corresponding matrix of anomalous dimensions is defined as 

A [n,p\[n',p'\ = ^[n,p\[n",p"]Fp.^[n",p"][n',p']- (5-3) 

A simple analysis of the diagrams shows that the matrix element -Z)n,p] [n'.p'j is proportional 
to h n ~ n , so that the elements with n < n! vanish (the parameter h = |h| appears only in 
the numerators of the diagrams; see Sec. 0. The elements with n = n' are independent 
of h and therefore they can be calculated directly in the isotropic model with h = 0. The 
block Z[ niP ] [ n yj can be then diagonalized by the changing to irreducible operators (scalars, 
vectors, and traceless tensors); but for our purposes it is sufficient to note that the elements 
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Z[ njJ ] [ n)P /] vanish for p < p' [the irreducible tensor of the rank p consists of the monomials 
with p' < p only, and therefore only these monomials can admix to the monomial of the 
rank p in renormalization]. Therefore, the renormalization matrix in Eq. ( |5.2| ) is triangular, 
and so is the matrix (|5.3|) . The isotropy is violated for h 0, so that the irreducible 
tensors with different numbers of the fields 9 can mix with each other even though their 
ranks are also different. In particular, the vector 8$ admixes to the irreducible tensor 
diQdjO — 5ij(d s 6d s 9)/d in the form of the traceless combination 2 5ij{h s 8 s 9)/d — hidj9 — hjdi9. 
In the following, we shall not be interested in the precise form of the basis operators, i.e., 
those having definite anomalous dimensions; we shall rather be interested in the anomalous 
dimensions themselves. The latter are given by the eigenvalues 7 [n, p] of the matrix ( |5.3|) , and 
in our case they are completely determined by the diagonal elements of the renormalization 
matrix ( ^2[ ): y[n,p] = V fl In Z [n , p] [n , p] . 

Now let us turn to the one-loop calculation of the constant ( |5.12| ) in the MS scheme. Let 
T(x; 6) be the generating functional of the 1-irreducible Green functions with one composite 


operator F[n,p] from Eq. (|5.1|) and any number of fields 6. Here x = (t, x) is the argument 
of the operator and 9(x) is the functional argument, the “classical analog” of the random 
field 9. We are interested in the 9 n term of the expansion of T(x; 9) in 9(x), which we denote 
r„.(x;d); it has the form 

T n (x]9) = ^J dxf-J dx n 9(x 1 )---9(x n )(F[n,p](x)9(x 1 )---9(x n ))i_ ir . (5.4) 

is represented diagramatically in the fol- 


In the one-loop approximation the function 
lowing manner: 


= F[n,p] + i 


(5.5) 


The first term is the “tree” approximation, and the black circle with two attached lines in the 
diagram denotes the variational derivative S 2 F[n,p]/8989. In the momentum representation 
it has the form 

d 2 F \n h 1 

r ( k ,= 89{k)89{q) = ~ lp (' p - X ) (<M • ■ • <V) {diOdiff) 1 


-Apl k h q s ( d i2 9 ■ ■ ■ d ip 9) d s 9 (.8$ 8$) l 1 
—2Z(k • q) ( d h 9---d lp 9) ( 8M9 ) z " 1 
-Al{l - 1 ) kfridjOdsO) (d h 9 ■ • ■ d lp 9) (8M9) 1 - 2 . 


(5.6) 


Strictly speaking, we had to symmetrize the right-hand side of Eq. (5.6) with respect to the 


indices i\ ■ ■ ■ i p and the momenta k, q. However, the symmetry is restored automatically 
after the vertex T(k, q) has been inserted into the diagram, which is why only one term of 
each type is displayed in Eq. ( |5.6|) and the required symmetry coefficients are introduced. 

The vertex (|5.6| ) contains (n — 2) factors of 89. Two remaining “tails” 9 are attached to 
the vertices 9'{wd)9 of the diagram (|5.5| ). It follows from the explicit form of the vertices 
that these two fields 9 are isolated from the diagram in the form of the overall factor 8989 ; 
cf. Sec. pd]. In other words, two external momenta, corresponding to these fields 9 , occur as 


an overall factor in the diagram, and the UV divergence of the latter is logarithmic rather 
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than quadratic; cf. the expression Q2.10|) , ( |3.11|) . Therefore, we can set all the external 
momenta and the “mass” m equal to zero in the integrand; the IR regularization is provided 
by the cut-off of the integral at q ~ m. Then the UV divergent part of the one-loop diagram 
can be written in the form 


dpOdid I — 


duj 

27T 


dq_ T (q ^ P p z(q) D v (u,q) 


(2n) d lu 2 + z/ 2 g 4 

The expression (|5.7|) is a linear combination of the integrals 


(5.7) 


Tij, p i — 


du 

27T 


dq QiQj P p i(q) D v (ui, q) 


( 2n) d lu 2 + u 2 q 4 

We perform the integration over t o and make use of the isotropy, namely, 


(5.8) 


d<i f{q) m 3 Ppi(q) = 


(d + 1 ) dpi djj SpiSij SpjSu 

d(d + 2 ) 


dqf(q) q 2 


This gives 


Tij, P i 


(d + 1 ) dpi (5 t j SpiSij SpjSu 

2 ud(d + 2 ) 


g^J 


(5.9) 


with the integral J from Eq. Q3.12|) . 

Substituting Eqs. (|5.6|) and (|5.9|) into Eq. ( |5.7|) gives the desired expression for the 
divergent part of the diagram (|5.5|) . In this expression we have to take into account all the 
terms proportional to the operator F[n,p] and neglect all the other terms, namely, the terms 
containing the factors of S tl i 2 etc. The latter determine non-diagonal elements of the matrix 
(m2), which we are not interested in here. Finally we obtain 


T„ ~ F[n,p] 


gp e JQ[n,p\ 
4 ud(d + 2) 


+ 


(5.10) 


where we have written 


Q np = 2n (n — 1) — (d + 1) (n - p) (d + n + p - 2) = 
= 2p {p — 1 ) — (d — 1 ) (n — p) (d + n + p). 


(5.11) 


The dots in Eq. (|5.10|) stand for the 0(g 2 ) terms and the structures different from F[n,p\, 
~ denotes the equality up to the UV finite parts; we also recall that n = p + 21. 

The constant 2y n ,p],[n,p] is found from the requirement that the renormalized analog T^ en = 
Z[“pi [ n p |T n of the function ( [5.101 ) be UV finite (mind the minus sign in the exponent); along 
with the representation (|3.13|) for the integral J and the MS scheme this gives the following 
result: 


Z 


[n,p] [n,p\ 


ag Q[n,p\ ^ (~u) s 
u 2(d - l)(d + 2 ) ^ e + si 7 ’ 


(5.12) 


with the polynomial Q[n,p\ from Eq. ( |5.11|) and the constant a is given in Eq. (|3 .1 5|) . 
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(5.13) 


For the anomalous dimension 


7 [n,p\ = 


we then obtain: 

agQ[n,p\ _ 

2u(u + l)(d-l)(d + 2) 


cf. Sec. |n| for the dimension 7 u . The critical dimension associated with the operator F[n,p] 
has the form A [n,p\ = 7 *[n,p]; see Eq. ( |4.19| ) and Table 1 ( 7 * denotes the value of 7 at the 
fixed point in question). For the nontrivial fixed points discussed in Sec. 13 we then obtain 

( = £ — r) for FPII, 

A [n,p\ = o/ - AAAA1 - ^ x { e/2 for FPIV, (5.14) 

r) = e/ 2 for FPV 


with the corrections of order 0(e 2 ). 

The expression ( |5.14| ) illustrates the general fact that the critical dimensions in the rapid- 
change and frozen regimes depend only on the exponents ( and e, respectively. It turns out 
that the dimension A [n,p\ at the point FPV is universal, i.e., it is independent of the free 
parameter w*, or, equivalently, of the specific choice of a fixed point on the curve described 
by Eq. (EII)- This is a consequence of the explicit form of the RG functions in the one-loop 
approximation (the same combination g/u(u + 1 ) enters into the beta functions and the 
anomalous dimension of the operator F[n,p]). We then expect that the exact dimension 
A [n,p\ at the point FPV is nonuniversal, and the dependence on w* will appear at the two- 
loop level. Another artifact of the one-loop approximation is the continuity of the dimension 
A[n,p\ at the crossover line e = 2r/ as a function of the exponents s, rj. 

The first-order result (|5.14|) for the operator F[ 2, 0] (the local dissipation rate) is in 
fact exact. The proof is based on certain Schwinger equation; it is almost identical to the 
analogous proof for the Kraichnan model, given in f32|, and will not be discussed here. 


The above analysis applies also to the case of a nonsolenoidal velocity field (compress¬ 
ible fluid). The transversal projector in Eq. (p~73|) is then replaced with k) + aQij(k), 
where Qij{ k) = kikj/k 2 is the longitudinal projector and a > 0 is an additional arbitrary 


parameter, the degree of compressibility. For the rapid-change regime ( |1.9|) , the dimension 
A[n,p] takes on the form 


A [n,p\ = 


< 


(d + 2 ) 


(n — p)(d + n + p) p(p 


1 ) (cx — 1 ) + a(n — p){n + p — 2 ) 


(d — 1 T ex) 


+ o( C 2 ), 

(5.15) 


in agreement with the p — 0 results obtained in Refs. [[H|] for the ‘tracer’ model and earlier 
in Ref. [^5|] for d = 1. In general case fll.7|) , additional superficial UV divergence emerges 
in the 1-irreducible Green function (9'6v) i_i r , and the second independent renormalization 
constant should be introduced as a coefficient in front of the new counterterm Q'(yd)Q. This 
case requires special analysis and will be discussed elsewhere^] (in particular, the nontrivial 
fixed point becomes infinite for the purely potential frozen velocity field, cf. |48],[49| for the 
random walks in random environment). 


2 See: N. V. Antonov, Anomalous scaling of a passive scalar advected by the synthetic compressible 
flow, |chao-dyn/990701S . 
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VI. OPERATOR PRODUCT EXPANSION AND THE ANOMALOUS SCALING 
FOR THE STRUCTURE FUNCTIONS AND OTHER CORRELATORS 


The representation (|4.22|) for any scaling function ((mr) describes the behavior of the 
Green f un ction for Ar >> 1 and any fixed value of mr. The inertial-convective range 
corresponds to the additional condition that mr « 1. The form of the function ((mr) is 
not determined by the RG equations themselves; in the theory of critical phenomena, its 
behavior for mr —> 0 is studied using the well-known Wilson operator product expansion 
(OPE); see, e.g., |45|j46|j 5l|| . This technique is also applicable to the theory of turbulence; 


see 


According to the OPE, the equal-time product Fi(x)F 2 (x') of two renormalized operators 
at x = (x + x')/2 = const and r = x — x' —> 0 has the representation 


F 1 (x)F 2 (x') = 


( 6 . 1 ) 


where the functions C a are the Wilson coefficients regular in m 2 and F a are all possible 
renormalized local composite operators allowed by symmetry, with definite critical dimen¬ 
sions A a . The renormalized correlator (Fi(x)F 2 (x')) is obtained by averaging Eq. ( 0 ) 
with the weight exp S ren , the quantities ( F a ) appear on the right-hand side. Their asymp¬ 
totic behavior for m —> 0 is found from the corresponding RG equations and has the form 
(F a ) (x m A “. From the operator product expansion m we therefore fold the following 


expression for the scaling function £(mr) in the representation (|4.22|) for the correlator 
(Fi(x)F 2 (x')): 


£( mr ) = y2,A a ( mr) 1 


( 6 . 2 ) 


where the coefficients A a = A a (mr ) are regular in (mr) 2 . 

In the models of critical phenomena, the leading contribution to the representations 
like (|6.2| ) is related to the simplest operator F — 1 with the minimal dimension A a = 0, 
while the other operators determine only the corrections that vanish for mr —> 0. The 
feature characteristic of the turbulence models is the existence of the so-called “dangerous” 
composite operators with negative critical dimensions [[32| . fl3|[39l |42|| . Their contributions to 
the operator product expansions determine the IR behavior of the scaling functions and lead 
to their singular dependence on m for mr —> 0. 

If the spectrum of the dimensions A a for a given scaling function is bounded from below, 
the leading term of its behavior for mr —> 0 is simply given by the minimal dimension. 
This is the case of the rapid-change model (see p2|,|33i), and, as we shall see below, of our 
model ( |2.4| ). [The exception is provided by the non-rapid-changes regimes, if the values of 
the exponents e, 7] are large enough. It is discussed in the subsequent Section.] 

Consider for dehniteness the equal-time structure functions of the scalar held: 


S n (r) = ([6>(t,x) - 6>(t,x')] n ), r = |x-x'|. (6.3) 

For the functions ( |6.3| ), the representation of the form (|4.22| ) is valid with the dimensions 
cIq = 0 and do = A G = nA g = —n. In general, the operators entering into the operator 
product expansions are those which appear in the corresponding Taylor expansions, and also 
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all possible operators that admix to them in renormalization. The leading term of the Taylor 
expansion for the function ( |3.3|) is given by the n-th rank tensor F\n, n] from Eq. (|5.1|) . The 
decomposition of F[n,Ji\ in irreducible tensors gives rise to the dimensions A [n,p] with all 
possible values of p; the admixture of junior operators gives rise to all the dimensions A [k,p\ 
with k < n. Therefore, the asymptotic expression for the structure function has the form 


S n (r) ~ (hr) 1 


n k 

EE 

k=0P=Pk 


[C kp ( mr) 


A [k,p] 


+ ••• 


(6.4) 


Here and below p k denotes the minimal possible value of p for given k , i.e., p k = 0 for k even 
and pk = 1 for k odd; C kp are some numerical coefficients dependent on e, rj, d, and on the 
angle between the vectors h and r. 

Some remarks are in order. 

The dots in Eq. ( |6.4| ) stand for the contributions of the order (mr) 2+ °^ and higher, 
which arise from the senior operators, for example, d 2 6d 2 6 or v 2 . 

In the original Kraichnan model, only scalar operators give contributions to the repre¬ 
sentations like (|6.4j) , because the mean values (F a ) of all the other irreducible tensors vanish 
owing to the isotropy; see f32| , |33|l . In the model (|2~4|), the traceless irreducible tensors ac¬ 
quire nonzero mean values, and their dimensions appear on the right-hand side of Eq. 0- 
In particular, the mean value of the operator 8^8j9 — 5ij(8 s 98 s 9)/d is proportional to the 
traceless tensor of the form 5ij(h s h s ) / d — hihj, its tensor indices are contracted with the 
indices of the corresponding coefficient C a in Eq. 0- 

The operators F[k,p] with k > n (whose contributions would be more important) do not 
appear in Eq. (|6.4j), because they do not appear in the Taylor expansion of the function S n 
and do not admix in renormalization to the terms of the Taylor expansion. 

The leading term of the expression (|6.4|) for mr —> 0 is obviously given by the contribution 
with the minimal possible dimension. The straightforward analysis of the explicit one-loop 
expression (|5.14j) shows that for fixed n, any d > 2, and any nontrivial fixed point, the 
dimension A [n,p\ decreases monotonically with p and reaches its minimum for the minimal 
possible value of p = p n . i.e., p — 0 if n is even and p — 1 if n is odd. Furthermore, this 
minimal value A [n,p n \ decreases monotonically as n increases, i.e., 


A[2 k, 0] > A[2 k + 1,1] > A[2 k + 2,0], 


[A similar hierarchy has been established recently in Ref. [^3| for the magnetic field advected 
passively by the rapid-change velocity in the presence of large-scale anisotropy.] Therefore, 
the desired leading term for the even (odd) structure function S n is determined by the scalar 
(vector) composite operator consisting of n factors 89 and has the form 


S n (r) oc ( hr) n (mr) A ^ n ’ Pri 


(6.5) 


with the dimension A [n,p\ given in Eq. (|5.14j) . 

For the rapid-change fixed point and even values of n, the total power of r in Eq. ( p.5| ) 
coincides with the exponent in the original isotropic Kraichnan model, calculated to the order 
0(() in [T7| and 0(1/d) in JTH] within the zero-mode approach, and to the order 0(0) in 
||32|| within the framework of the RG. We also note that the anomalous dimensions associated 
with the operators F[2k,2] were calculated in |32j to the order 0(C 2 ); the exact dimension 
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of the operator F[ 2, 2] was found in |14|. [It should be noted that the decomposition of the 
total exponent in Eq. (|6.5|) into the critical dimension of the composite operator and the 
critical dimension of the structure function itself differs from the analogous decomposition 
for Kraichnan’s model, as a result of the difference in canonical dimensions; see Sec. 0 . 

The result ( |S . 5| ) for the third-order structure function in the rapid-change model coincides 
with the 0(C) result obtained in | 2 ^| within the zero-mode technique; see also the earlier 
paper |T3 for the three-dimensional result. We note that the exponents — 7C/5 and 3C/5 


from [0 should be identified with the anomalous dimensions A[3,1] and A[3, 3], respectively. 
The result (To) for n = 3 is also in agreement with the 0(l/d) result obtained in (jHJ, with 
the identification 7 + 1 — A = 3 + A[3,1]. 

For the case of the frozen velocity field (FPIV), the first-order results for the even struc¬ 
ture functions were presented in j32|. We also note that they satisfy the exact inequalities 
obtained for the time-independent case in [|38| . 

The analysis given above is extended directly to the case of other correlation functions. 
For example, the analog of the expression (|6.5|) for the equal-time pair correlation function 
of the operators (|5.1|) has the form 


(F[n,p] F[n',p ']> ~ h n+n ' (A r )- A M-A[n' )Pn ,]( mr )A[n+n',p n+n ,] 


( 6 . 6 ) 


Some special cases of the relation (| 6 . 6 |) for the rapid-change model were obtained earlier in 

Refs. 0 00 

Another interesting example is the equal-time pair correlator (9 n {t 1 x.)9 k {t,x.')). Substi¬ 
tuting the relations = 0 and do = A G = — (n + k ) into the general expression (|4.22|) 
gives (9 n 9 k ) = r n+k ^(rnr), and the small mr behavior of the scaling function £(mr) is found 
from Eq. ( |6.2| ) (here and below, we do not display the obvious dependence on h ). In 
contrast to the previous example, composite operators in the expansion (ED can involve 
the field 9 without derivatives. The leading term in Eq. ( |6.2| ) is then given simply by 
the operator 9 n+k with Ap = —(n + k), while the first correction is related to the mono¬ 
mial ( di9di9)9 n+k ~ 2 whose critical dimension is easily found to be Ap = — (n + k) + A u 
with A w from Eq. (|4.20|) . Therefore, in the inertial range our correlator has the form 


( 9 n 9 k ) ~ Ci'mT^+O — C 2 m"^ + ^(mr) A “ + large constant minus a powerlike correction 
(the signs of the constants <7 are explained by the fact that the correlator is positive and 
decreases as r grows). In the structure functions ( |Q|) all the contributions related to oper¬ 
ators containing fields without derivatives cancel out to give the behavior ( |6.4| ) , determined 
by the operators constructed only of held derivatives. 

Finally, we note that the hierarchy of critical dimensions A [n,p\, established in Sec.[V]. 
persists also for the nonsolenoidal velocity held, see (|5.15|) . Therefore, the asymptotic ex¬ 
pressions like 0 , 0 and (| 6 . 6 | ) remain valid for the compressible case (a ‘tracer’ in 
terminology of 0 ) with the exponents A [n,p] given in (|5.15|). 


VII. SUMMATION OF THE DANGEROUS CONTRIBUTIONS FROM THE 
POWERS OF THE VELOCITY FIELD 


A new interesting problem emerges as the parameters e and 7 increase and the velocity 
held becomes dangerous; see Eq. ([4.21|). Owing to the Gaussianity, all its powers also 
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become dangerous with the dimensions A[n ?1 • • • v ln ] = nA v . The analysis of the diagrams 
shows that in the rapid-change regime, these operators do not contribute to the operator 
product expansions of the equal-time correlators like ( |6.3| ) or (|A6|), but all the contributions 
of the scalars (v 2 ) n do appear in those OPE for the non-rapid-change regimes. The spectrum 
of their dimensions is unbounded from below, and in order to End the small mr behavior 
we have to sum up all their contributions oc (rnr) nA ^ in the representation ( |6.2|) . We have 
employed the infrared perturbation theory in the form developed in to perform the 

required summation for the structure function S 2 in the frozen regime and within the one- 
loop approximation for the Wilson coefficients. In this case, the velocity becomes dangerous 
for e > 2 (region FPIV& on Fig. I). 

The function S 2 is represented in the form 


S 2 = J Vi [e(t, x) - Hit x')] 2 exp S(4>) 


(7.1) 


with the action functional from Eq. (|2.4| ). Following we split the velocity held in Eq. 

( ffTTp into two components, v(ic) = v < (x) + v > (x), referring to the “soft” component, v<, 
all the Fourier modes with k < fc*, and to to the “hard” component, v>, all the remaining 
modes with k > k*. Here fc* is a fixed arbitrary separating scale, which will not enter into 
the final expressions. Since we are interested in only the contributions of the operators 
(v 2 ) n into the OPE, we can neglect the spacetime inhomogeneity of the soft held. It then 
becomes a random variable (rather than a random held) with the statistics determined by 
the relation 


(v< • •■v<) = (v(x) • • • v(x)). 


(7.2) 


Furthermore, we confine ourselves to the one-loop approximation for the corresponding 
Wilson coefficients, so that we can omit the contribution of the hard held in the vertex 
0'(vd)0. Then the integration over the helds 6, O ', and v> in Eq. (|7.1|) gives: 


S 2 = 2 


[■ dk 

r duj 

1 (2v r) d J 

(2x)[ 


1 — exp(ik • r) 


P(k)A,(u;,k) 


u; 


k ) 2 + z/q/c 4 ’ 


(7.3) 


where P(k) = hihjPij( k), the correlator D v is given by Eq. 0, and the averaging over v< 
with the statistics ( |7.2|) is to be performed. The mean values (|74^) in Eq. (|7.3|) correspond 
to the contributions from (v(x) ■ ■ ■ v(x)) in the in the representation (|6.2|) for S 2 . For the 


rapid-change model, the correlator D v is independent of the frequency; see Eq. (|1.9| ) . It then 
follows from the expression (|7.3|) that the dependence on v< vanishes after the integration 
over u>, which means that the operators (v 2 ) n give no contribution to the OPE for S' 2 . Now 
let us turn to the case of the frozen velocity held. We substitute the correlator (|1.10|) into 
Eq. O and perform the integration over a;; this gives: 


S2 = g'^o J 


dk 

(2^L 


1 — exp(ik • r) 


P(k) (k 2 + m 2 ) d / 2+l £ / 2 
(v< ■ k ) 2 + z/g/c 4 


(7.4) 


A brief digression is required here. We are interested in the small mr behavior of the 
scaling function ((mr) from Eq. ( |6.2|) , so that we have to combine the expression (|7.4| ) with 
the representation (|4.22| ) for the function S 2 . In renormalized variables, the latter is written 
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in the form S 2 = ( hr) 2 f (fir , g", u, mr ), where / is some function of completely dimensionless 
arguments. The function £(mr) is then given by the relation ((mr) = f(l,g”,0,mr) with 
g” from Eq. (|4.12| ) (we recall that for the frozen regime, w* = 0). The renormalization of 
Eq. (|7.4|) in our approximation reduces to the replacement g{[ —> g"ii £ , u 0 —> u, and the 
changeover to the scaling function is given by the substitution g" —* g", fi —> 1 jr. From 
now on, all these substitutions are implied. The expansion of the denominator in Eq. ( |7.4|) 
in (v< • k ) 2 gives: 


---= — y'(-ir— 

(v< • k ) 2 + v 2 k A u 2 k A j ~ ) 


k) 


2 n 


(7.5) 


v 2n k An 

It follows from Eqs. Q1.10| ) and (|7.2| ) that the correlators of the soft field have the form 

((u<)u ■'' («<)«») = D [S ili2 S i3i4 ■ ■ ■ S i2n _ li2n + all possible permutations], (7.6) 

where 

elk 


D ~ r £ v 2 


J (2 ) d d ^ 2+1 £ ^ 2 — v 2 m 2 (mr) £ . 


(7.7) 


(here and below ~ denotes the equality up to a numeric factor). Strictly speaking, the 
integral ( fT7| ) should be cut off from above at k ~ fc*. For £ > 1, the cut-off can be removed 
without changing the singularity at m -> 0. The averaging Eq. (El) over v< gives: 


E 


(2 n) 


D 


z/2/c 4 — 0 n\ 


l-Z) 


z = 


2 k 2 v 2 


m 2 (mr)~ 


(7.8) 


[we note that (2u)!/2 n n! = (2 n — 1)!! is the number of terms in Eq. (|7.6|) ]. The small mr 
limit implies z —> 00 . The large z behavior of the series in Eq. ( |7.8| ) is found from the 
following integral representation: 


E 

n =0 


(2n)! 


n\ 


(—z) n = / dt exp(—zt 2 — t) = J-k/Az l + 0(l/y/z) 


(7.9) 


Substituting Eq. ( |7 into Eq. ra gives: 


S 2 — rn 2 {mr) £//2 J dy 


1 — exp(iy • mr) 


P(y)y- 3 {l+y 2 )- d/2+1 - £/2 , 


(7.10) 


where we have changed to the dimensionless integration variable y defined so that k = my. 
Expansion of the quantity in the square brackets for small mr gives 


S 2 ~ (mr) - Vr . j d.yP(y)y i y j y 3 (1 + j/ 2 ) 1,2 


(7.11) 


[the first term of the expansion gives no contribution to Eq. (|7.11|) , owing to the evenness of 
the rest of the integrand]. The vector indices can be isolated from the integral ( 7.1 1|) ; this 
gives rise to the two structures, 5ijh 2 and /'q/'q, with the scalar coefficients proportional to 
the integral / dyy _1 (l + y 2 )~ d / 2+1 ~ e / 2 . One can easily check that for £ + g > 1, it is both 
IR and UV convergent, so that the leading terms of the small mr behavior of the integral 
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( [T.1C| ) are indeed obtained simply by the expansion of the integrand and have the form 
h 2 r 2 (mr)~ £ t 2 and (h • r ) 2 (mr) _e / 2 . 

Therefore, it turns out that the contributions of the operators (v 2 ) n sum up to the 
powerlike expression r 2 (mr) _e / 2 . In other words, the infinite sum of these dangerous oper¬ 
ators gives to the function S 2 contribution of the same form as the single operator F[2,0], 
and therefore the IR behavior of S 2 is given by the same expression (|6.5|) for all values of 
the exponent e [of course, the corresponding amplitudes for e > 2 acquire an additional 
contribution from the operators (v 2 ) n j. 

The infinite family of the dangerous operators (v 2 ) n also occurs in the RG approach 
to the stochastic NS equation; see |39|--|42|| . In that case, their summation in the OPE for 


different-time correlators leads to the expressions that are analogous to |47j and describe the 
well-known sweeping effects. The contributions of these operators vanish when one changes 
to the equal-time Galilean invariant functions, for example, the structure functions of the 
velocity field. In this connection, it should be stressed that the appearance of the operators 


r 2\n 


in the structure functions of the model ([2.4]) is not an artifact of the synthetic velocity 
statistics. One can check that in the presence of a mean gradient, the same effect occurs 
even though the velocity is generated by the Galilean covariant stochastic NS equation. On 
the contrary, if the effective scalar noise in the diffusion equation is proportional to the delta 
function in time (as in original Kraichnan’s model), the operators (v 2 ) n are absent from the 
equal-time correlators whatever be the velocity statistics. One can probably consider this 
effect as an additional source of the breakdown of the Kolmogorov-Obukhov theory for the 
passive scalar advection. 

The summation given above can be viewed as a possible model of the origin of the 
anomalous scaling in the structure functions of the velocity field: it was argued in [50 


that the singular dependence on mr of the equal-time correlators for the stochastic NS 
equation can be related to infinite families of dangerous operators. The summation can also 
be performed in a more formal way, without referring to the infrared perturbation theory, 
by using only the operator product expansion for the quantity [9(t, x) — 9(t, x ')] 2 and not 
only its average (ED- It is also possible to take into account all the composite operators 
constructed of the velocity and its time derivatives; see |54j . 


VIII. “EXOTIC” SCALING REGIMES 


So far, we have considered the standard formulation of the asymptotic problem, for which 


the dimensional bare charges entering into the velocity correlator (|1.7| ) depend on the UV 
scale only; see Eq. (|1.6|). I 11 a number of papers, e.g., [I].[T|[TI] [T| , a different version of 
the problem was considered, in which the velocity correlator depends also on the IR scale 
L = m~ 1 . The RG method is applicable also to this “exotic” (from the viewpoints of 
the theory of critical behavior) situation. No additional calculation of the fixed points or 
anomalous exponents is required; only the regions of stability of the fixed points FPI-FPV 
can change .0 


'^This Section is not included into the journal version. 
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As already said in Sec. [TV], solutions of the RG equation (|3.8|) can be represented in terms 
of invariant charges g , u. The relations ( 4.17 ) determine implicitly the invariant charges as 
functions of the bare parameters; the former depend on the scales A, m only through their 


dependence on the latter. The solution of the equations (|4. 1 7|) for any given g 0 , u 0 would 
show which of the fixed points is “chosen by the RG flow.” However, the complete solution 
of this problem is too difficult to achieve; the main obstacle is that the renormalization 
constants in Eq. ( |4.17p must be taken in the one-loop order of the RG, while we know them 
only in the one-loop order of the ordinary perturbation theory; see Eq. (|3.16|) . In other 
words, we know the renormalization constants only up to the terms of order 0(g), and here 
we need to take into account all the terms of the form (g/e) n etc. 

We shall not try to solve the problem completely, we rather give a simple recipe that 
allows one to distinguish between the nontrivial scaling regimes. The renormalization con¬ 
stants can be eliminated from Eq. ( [4.17|) using the relation (|3.5|) ; this gives 


x(k) = 4,= ^ 


u 




( 8 , 1 ) 


with ( = e — g, see Eq. ( |l.jj| ). Of course, the identity ( |3.1| ) contains less information than 
the full system ( |4.17|) ; in particular, we cannot judge from Eq. (|S.1|) whether or not the 
invariant charges in the IR regime are attracted by one of the nontrivial fixed points FPII, 
FPIV or FPV. Fortunately, if we know (or assume) that this is the case, we can definitely 
say which point is involved. Indeed, for FP11 we have u —> oo, g/u 2 = const (see Sec. 
B and therefore x(k) —> 0 at this fixed point. For FPIV (u —> 0, g/u = const) we have 
x(k) —> oo and for FPV the function x(k) remains 0(1). We also note that for this case 
x(k) is independent of k owing to the equality ( = g. Therefore, taking the limit k —> 0 in 
( jg.l| ) leads to the exact relation go/u^ = g*/u 3 for FPV, which was used in Sec. [TV| when 
deriving Eq. ( 4.1(j| ). 

For the standard regime ( |1.6P it follows from Eq. (|8.1|) that x(k) oc (k/ A)^. Then the 
above considerations show that in the IR asymptotic region, k « A, the invariant charges 
approach the point FPII for ( — g < 0, FPIV for ( — g > 0, and FPV for ( = g, in agreement 
with the analysis given in Sec. 0- 

Now let us turn to the example of the “exotic” regime for which the characteristic 
frequency of the velocity held depends only on the IR scale, u oc (Wk 2 ) 1 ^ 3 (k /m) 4 / 3_r? [the 
dependence on the mean energy dissipation rate, W, is implied], while the equal-time second- 
order structure function of the velocity is expressed only through the UV scale, 5*2(r) oc 
(Wr) 2 / 3 (Ar) -8 / 3+ G | " ?? . The comparison with Eq. 0 ) gives uqVq oc W 1 ^ 3 (k/m) 4 0 v and 
go^o/uo oc iy 2 / 3 A -8 / 3 "G +r? . From Eq. ( [S.l|) we then obtain 


X(k) oc (A/k)“ (k/mf, 


( 8 . 2 ) 


with the exponents 


a = -8/3 + C + g, f3 = - 8/3+ 2g. (8.3) 

It follows from Eqs. (|8.2|) , ( g.3P that the regions of stability of the fixed points have changed: 
the point FPII is approached if a, (3 < 0, FPIV is approached if a, (3 > 0, and the nonuni- 
versal regime FPV takes place for the unique choice a = (3 = 0 that exactly corresponds to 
the Kolmogorov velocity correlator (£ = g = 4/3). 
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The most interesting scaling regime arises when the signs of the exponents a and (3 are 
opposite, say, a > 0, /3 < 0. Then near the upper bound of the inertial range, k » m, 
A ~ fc, we have x(fc) << 1, and the critical dimensions are given by the rapid-change fixed 
point FPIV; near the lower bound of the inertial range, A >> k, k ~ m, we have y(/c) << 1, 
and the “frozen” point FPV works. Therefore, the same critical regime is described by two 
different fixed points and, consequently, two sets of critical dimensions! ft is tempting to 
relate this fact to the effect observed experimentally in the boundary layer: the second-order 
structure function exhibits two power-law regimes with a pronounced break in the exponent 
||59|| . Of course, one should not insist too much on this interpretation. 


IX. CONCLUSION 


We have applied the RG and OPE methods to the simple model describing the advection 
of a passive scalar by the synthetic velocity held and in the presence of an imposed linear 
mean gradient. The statistics of the velocity is Gaussian, with a given self-similar correlator 
with finite correlation time. 

We have shown that the model possesses the RG symmetry, and the corresponding RG 
equations have several fixed points. As a result, the correlation functions of the scalar held in 
the inertial-convective range exhibit various types of scaling behavior: diffusive-type regimes 
for which the advection can be treated within the ordinary perturbation theory, and three 
nontrivial convection-type regimes for which the correlation functions of the model reveal 
anomalous scaling behavior. The stability of the hxed points (and, therefore, the choice of 
the scaling regime) depends on the values of the two exponents e and rj, entering into the 
velocity correlator. 

The explicit asymptotic expressions for the structure functions and other correlation 
functions in any space dimension are obtained; the anomalous exponents are calculated 
to the first order of the corresponding £ expansions. For the first nontrivial regime the 
anomalous exponents are the same as in the rapid-change version of the model; for the 
second they are the same as in the model with time-independent (frozen) velocity field. In 
these regimes, the anomalous exponents are universal in the sense that they depend only 
on the exponents entering into the velocity correlator; what is more, they depend on the 
only exponent (( = £ — rj and e) remaining in the corresponding limit. For the last regime 
the exponents are nonuniversal: in principle, they depend also on the values of the coupling 
constants. It turns out, however, that they can reveal the nonuniversality only in the second 
order of the e expansion. 

A serious question is that of the validity of the e expansion and the possibility of the 
extrapolation of the results, obtained within the e expansions, to the finite values e = 0(1). 
In Refs. [^2| and [28|], the agreement between the nonperturbative results and the £ expansion 
has been established on the example of the triple correlation function in the rapid-change 
model. In particular, in [|28| the exponent A [3,1] (we use the notation introduced in Sec. 
0 has been calculated numerically for all 0<C<2 within the zero-mode approach. It was 
shown that for small (, this nonperturbative result agrees with the expansion in (, while for 
C = 2 it coincides with the exact analytic result A [3,1] = —2 obtained previously in [ED . 
In the paper |25|, the one-dimensional version of the rapid-change model has been studied 


both numerically and analytically within the zero-mode approach; the analytic expressions 


29 













for the anomalous exponents obtained within the £ expansion have also been found to agree 
with the nonperturbative numerical results. Finally, in Ref. [29] the analytic 0(e) result has 


been confirmed by the numerical experiment on the example of the fourth-order structure 
function in three dimensions. 

In this connection, we also note that a number of exact analytic results appear to be in 
agreement with the corresponding e expansions: the exponent A[2, 2], calculated exactly in 
0, the exponent for the second-order structure function of a passively advected magnetic 
held [^7], and the second-order exponent for a scalar advected by the nonsolenoidal (“com¬ 
pressible”) velocity held ]33[; the corresponding expansions in £ (to the order ( 2 ) have been 
calculated within the RG approach in |32|J33| . These facts support strongly the applicability 


of the e expansion, at least for low-order correlation functions. 

In the paper |ll], a closure-type approximation for the rapid-change model, the so-called 
linear ansatz, was used to derive simple explicit expression for the anomalous exponents for 
any 0 < ( < 2, d, and n, the order of the structure function. For £ = 1, the predictions of the 
linear ansatz appear to be consistent with the numerical simulations |T^j23|j29|j30|l ; they are 
also in agreement with some exact relations fT3| , |2T| , |23| . However, they do not agree with the 
results obtained within the zero-mode and RG approaches in the ranges of small (, 2 — ( or 
1/d. In fact, there is no formal contradiction between the perturbative results and the linear 
ansatz: the violation of the latter in the aforementioned limits can be related to the fact 
that they have strongly nonlocal dynamics in the momentum space; see, e.g., the discussion 
in Refs. l[Ri| . [2b|l . On the other hand, the numerical divergence of the predictions given by the 
linear ansatz and £ expansion for the fourth-order structure function at ( ~ 1 is, roughly 
speaking, of the same order of magnitude as the difference between the nonperturbative 
numerical results and the perturbative small-C results for the triple correlator, as one can see 
from the figures presented in One can think that the series in obtained within the RG 
or zero-mode approaches, give correct formal expansions of the (unknown) exact exponents, 
while the linear ansatz gives a good approximate expression for the same quantities near 
( = 1. We also note that the numerical agreement between the expansion in ( and the exact 
results is expected to worsen as n increases, because the actual expansion parameter is n( 
rather than ( itself; see |32| , |33| . 

Of course, it is not impossible that new dangerous operators arise for some finite value 
of the RG expansion parameter. The example is provided by the frozen regime of the model 
in question (see Sec. jVTl j). In principle, this effect can lead to a qualitative changeover in 
the IR behavior of the correlators as £ increases, but in our model this is not the case, at 
least for the second-order structure function. 

We also note that the first-order expressions for the anomalous exponents in Eq. (|5.14 ) 
look alike for the rapid-change and frozen regimes, but we expect that the analytic properties 
of the series are different. The aforementioned exact results suggest that in the rapid-change 
models, the series in ( have finite radii of convergence. In the language of the held theory, 
this is related to the fact that in the rapid-change models, there is no factorial growth of the 
number of diagrams in higher orders of the perturbation theory: a great deal of diagrams 
vanishes owing to retardation and the fact that the velocity correlator contains the delta 
function in time. This is no longer so for the regimes in which the correlation time remains 
finite, and we expect that the series in £ for these regimes are only asymptotical, as in most 
models of the critical behavior. 
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Let us conclude with a brief discussion of the passive advection in the non-Gaussian 
velocity held governed by the nonlinear stochastic NS equation. In this case, one has to 
add the nonlinear term (vjdj)vi to the left-hand side of the equation (|1.2| ) and set 7] = 0 
in Eq. (|1.4j) . The RG approach is also applicable to this model; the analysis of the UV 
divergences shows that the basic RG functions are the same as for the model with h = 0. 
The RG analysis of the latter has been accomplished in |h|. It shows that the model 
possesses a nontrivial IR stable fixed point; its coordinate has been calculated in [B3J in the 
first order of the e expansion .[] The inclusion of a nonzero mean gradient h ^ 0 gives rise 
to anomalous scaling; the analysis given in Sec. 0 can also be extended to this case. For 
small e, the anomalous exponents are given by the relation (|5.13| ), in which one should take 
g/u(u + 1) = e/3 a at the fixed point, with the coefficient a from Eq. (|3.15|) (we use the 
notation introduced above; the definition of the parameters e, a, and u in [Q is slightly 
different). Despite the non-Gaussianity, the critical dimensions of the powers of the velocity 
held are given by the simple linear relation A [v il ■ ■ ■ v ln ] = nA v = n(l — e/3); see 
(in the notation of the papers 


0 


39|-f42|, e should be replaced with 2e). Therefore, all these 


operators are dangerous for £ > 3, and the summation of their contributions is required. For 
the different-time correlators, it has been accomplished in ||39|.[40|| ; for the structure functions 


it can be performed in the one-loop approximation as in Sec. IV III above and leads to an 


analogous conclusion: the behavior of the second-order structure function does not change 
for £ > 3. For £ > 4, the composite operator of the local energy dissipation rate also becomes 
dangerous |52|, possibly along with all of its powers [|57|]; some other dangerous operators 
arise for e > 6 and further |54|,|58l|. The identification of all the other dangerous operators 


and summation of their contributions in the operator product expansions remains an open 
problem. 


ACKNOWLEDGMENTS 

I have benefited from discussions with L. Ts. Adzhemyan, M. Hnatic, J. Honkonen, 
M. Yu. Nalimov and A. N. Vasil’ev. I am thankful to V. S. L’vov, A. Mazzino and M. Ver- 
gassola for useful comments on the paper subject. 

The work was supported by the Russian Foundation for Fundamental Research (Grant 
No. 99-02-16783) and by the Grant Center for Natural Sciences of the Russian State Com¬ 
mittee for Higher Education (Grant No. 97-0-14.1-30). 


! The results obtained in 55] were also rederived later in Refs. [561. 
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Fig.I. Regions of stability for the fixed points in the model (2.4). 
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TABLES 


TABLE I. Canonical dimensions of the fields and parameters in the model (p.4[). 


F 

9 

9' 

V 

v, v 0 

m, M, fi, A 

9o 

u 0 

g , u, h 

d k F 

-1 

d+ 1 

-1 

-2 

1 

e + g 

g 

0 

d F 

0 

0 

1 

1 

0 

0 

0 

0 

dp 

-1 

d + 1 

1 

0 

1 

e + g 

g 
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